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We study the lepton forward-backward asymmetry Afb and the longitudinal K* ** polarization Fl, 
as well as an observable P 2 derived from them, in the rare decays B —>■ K*£'^£~, where £'^£~ is either 
e''"e“ or pAp.~, using the full sample of 471 million BB events collected at the T(4S') resonance 
with the BABAR detector at the PEP-II e^e“ collider. We separately fit and report results for the 
K*^(892)£~^£~ and K*~^ (892)£~^£~ final states, as well as their combination K*£~^£~, in five disjoint 
dilepton mass-squared bins. An angular analysis of B'^ —>■ K*'^£'^£~ decays is presented here for 
the first time. 

PACS numbers: 13.20.He, 12.15.-y, ll.30.Er 


I. INTRODUCTION 


The decays B -A K*{892)i'^£~, where K* -A Kn 
(hereinafter, unless explicitly stated otherwise, K* refers 
generically to the Ar*(892)) and £~^£~ is either an e+e“ 
or pA pair, arise from flavor-changing neutral-current 
(FCNG) processes, which are forbidden at tree level in 
the Standard Model (SM). The lowest-order SM pro¬ 
cesses contributing to these decays are the photon pen¬ 
guin, the Z penguin and the W'^W~ box diagrams shown 
in Fig. [T] Their amplitudes are expressed in terms of 
hadronic form factors and perturbatively-calculable ef¬ 
fective Wilson coefficients, Cf®, C|® and G®®, which 
represent the electromagnetic penguin diagram, and the 
vector part and the axial-vector part of the linear combi¬ 
nation of the Z penguin and W~^W~ box diagrams, re¬ 
spectively [ll-0|- Non-SM physics may add new penguin 
and/or box diagrams, as well as possible contributions 
from new scalar, pseudoscalar, and/or tensor currents, 
which can contribute at the same order as the SM dia¬ 
grams, modifying the effective Wilson coefficients from 
their SM expectations An example of a non-SM 

physics loop process is shown in Fig. O other possible 
processes could involve e.g., non-SM Higgs, charginos, 
gauginos, neutralinos and/or squarks. As a function of 
dilepton mass-squared = rrij +^_, the angular distri¬ 
butions in H —>■ K*Ai~ decays are notably sensitive to 
many possible sources of new physics, with several collab¬ 
orations presenting results over the past few years [11- 



FIG. 1: Lowest-order SM Feynman diagrams for h s£'^£ . 


At any particular value, the kinematic distribution 
of the decay products of B —5> K*Ai~ and the CP- 
conjugate B -A K*Ai~ process depends on six transver- 
sity amplitudes which, neglecting GP-violating effects 
and terms of order and higher, can be expressed as 
a triply differential cross-section in three angles: 9k, the 
angle between the K and the B directions in the K* 
rest frame; 9i, the angle between the A{£~) and the 
B{B) direction in the Ai~ rest frame; and (p, the angle 
between the A£~ and Kir decay planes in the B rest 
frame. From the distribution of the angle 9k obtained 
after integrating over (p and 9i, we determine the K* lon¬ 
gitudinal polarization fraction using a fit to cos 9k of 
the form 


*Now at: University of Tabuk, Tabuk 71491, Saudi Arabia 
^Now at: Laboratoire de Physique Nucleaire et de Hautes Energies, 
IN2P3/CNRS, F-75252 Paris, France 

^Now at: University of Huddersfield, Huddersfield HDl SDH, UK 
^ Deceased 

^Now at: University of South Alabama, Mobile, Alabama 36688, 
USA 

**Also at: Universita di Sassari, 1-07100 Sassari, Italy 


1 dr 

r{q^) d(cos6»*-) 


^FL{q'^) cos^ 9 k + 

^(1 - FL{q^)){l - cos^ 9k) ■ (1) 


We similarly determine the lepton forward-backward 
asymmetry Afb from the distribution of the angle 9i 
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FIG. 2: Feynman diagram of a non-SM Higgs penguin pro¬ 
cess. 


obtained after integrating over (j) and 9k, @ 

r(kd(^ = + 

^(l-ei(9"))(l+':os=»,) + 

ApBiq"^) COS9i. ( 2 ) 


We ignore here possible contributions from non-resonant 
S-wave B —>■ Ktt£^£~ events. The rate for such events 
has been shown to be consistent with zero [1^, with an 
upper limit (68% CL) across the entire dilepton mass- 
squared range of < 4% of the B K*{Ktt)£~^£~ branch¬ 
ing fraction [^ . The presence of an S-wave component 
at this level was shown to lead to a relatively small abso¬ 
lute bias on the order of 0.01 for and Afb', this small 
bias is ignored here given the relatively larger magnitude 
of our statistical and systematic uncertainties. Essen¬ 
tially no contributions from low-mass tails of the higher 
K* resonances are expected in the iL*(892) mass region 
considered here. 

We ignore small g^-dependent theory corrections in 
the large-recoil ^ 2 GeV^/ A region given the current 
experimental uncertainties on the angular observables, 
which are relatively large compared to these small cor¬ 
rections in the underlying SM theory expectations Q. 
We determine Fl and Afb va the five disjoint bins of 
defined in Table HI We also present results in a q^ range 
1.0 < ( 7 q < 6.0GeV^/c'*, the perturbative window away 
from the q^ ^ 0 photon pole and the cc resonances at 
higher where theory uncertainties are considered to 
be under good control. An angular analysis of the de¬ 
cays B'^ K*^£'^£~ is presented here for the first time. 

We additionally present results for an observable derived 
from Fl and Afb, P 2 = (— 2/3)*Afs/(1 —Aj;,), with less 
theory uncertainty, and hence greater sensitivity to non- 
SM contributions, than either Fl or Afb alone [^.1^. 


II. EVENT SELECTION 

We use a data sample of ^ 471 million BB pairs, corre¬ 
sponding to 424.2 ± 1.8fb“^ [13, collected at the r(4S') 


resonance with the BaBAR detector [U at the PEP-II 
asymmetric-energy e+e“ collider at the SLAG National 
Accelerator Laboratory. Charged particle tracking is pro¬ 
vided by a five-layer silicon vertex tracker and a 40-layer 
drift chamber in a 1.5 T solenoidal magnetic field. We 
identify electrons and photons with a CsI(Tl) electro¬ 
magnetic calorimeter, and muons using an instrumented 
magnetic flux return. We identify charged kaons using a 
detector of internally reflected Cherenkov light, as well 
as AE/dx information from the drift chamber. Charged 
tracks other than identified e, /i and K candidates are 
treated as pions. 

We reconstruct B -A K*£~^£~ signal events in the 
following final states (charge conjugation is implied 
throughout unless explicitly noted): 

• B^ -A K*^{^ KgTT^)fi^fj,~; 

• B° ^ K*°{-a K+Tr-)fi+n-] 

• B^ -A Ar*+(—>■ K^TT^)e'^e~; 

• B^ -A A'*+(— KgTT^)e'^e~; 

• B° -)■ K*°{-a K+F-)e+e-. 

We do not include the decays —>■ A'*+(—>■ 

Ar+7r°)/x+/i“ and B^ -A K*^{-a KgTT^)£~^£~ in our anal¬ 
ysis. The expected signal-to-background ratio for these 
final states relative to the five chosen signal modes listed 
above is very poor, with ensembles of pseudo-experiments 
showing that inclusion of these extra modes would yield 
no additional sensitivity. 

We require K* candidates to have an invariant mass 
0.72 < m{KTT) < l.lOGeV/c^. Electron and muon can¬ 
didates are required to have momenta p > 0.3 GeV/c in 
the laboratory frame. The muon and electron misiden- 
tification rates determined from high-purity data control 
samples are, respectively, ~ 2% and < 0.1% [3l|, and 
backgrounds from particle misidentification are thus sig¬ 
nificant for B -A K* pA pL~ candidates only. We combine 
up to three photons with an electron candidate when 
the photons are consistent with bremsstrahlung from the 
electron. We do not use electrons that are associated 


TABLE I: Definition of the bins used in the analysis. The 
nominal B and K* invariant masses are given by ms and 
rriK*, respectively. 


bin 

q^ min (GeV^/c'^) 

q^ max (GeV^/c'*) 

91 

0.10 

2.00 

92 

2.00 

4.30 

93 

4.30 

8.12 

94 

10.11 

12.89 

95 

14.21 

(ruB - mK*A 

9o 

1.00 

6.00 
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with photon conversions to low-mass e+e“ pairs. We 
reconstruct candidates in the tt+tt” final state, re¬ 
quiring an invariant mass consistent with the nominal 

mass, and a flight distance from the e+e“ interaction 
point that is more than three times the flight distance 
uncertainty. Neutral pion candidates are formed from 
two photons with > 50 MeV, and an invariant mass 
between 115 and 155MeV/c^. In each final state, we uti¬ 
lize the kinematic variables rriES = ~ P*b 

AE = Eg — EcM/2, where p*g and Eg are the B momen¬ 
tum and energy in the T(45') center-of-mass (CM) frame, 
and Ecm is the total CM energy. We reject events with 
Toes < 5.2GeV/c^. 

To characterize backgrounds from hadrons misidenti- 
fied as muons, we study candidates, where h is a 

charged track with no particle identification requirement 
applied. We additionally use a K*e^pL^ sample, where 
no signal is expected because of lepton-flavor conserva¬ 
tion, to model the combinatorial background from two 
random leptons. For both e+e“ and p.~ modes, we 
veto the J/if) (2.85 < mi+i- < 3.18GeV/c^) and 
(3.59 < mg+£- < 3.77 GeV/c^) mass regions. These ve¬ 
toed events provide high-statistics control samples of de¬ 
cays to final states identical to the signal modes here that 
we use to validate our fitting procedures. 

Random combinations of leptons from semileptonic 
B and D decays are the predominant source of back¬ 
grounds. These combinatorial backgrounds occur in both 
BB events (“BR backgrounds”) and e+e“ —?> qq con¬ 
tinuum events (“gq backgrounds”, where q = u,d,s,c), 
and are suppressed using eight bagged decision trees 
(BDTs) [13 trained for suppression of: 

• BB backgrounds in e+e“ modes at low q"^; 

• BB backgrounds in e+e“ modes at high 9 ^; 

• BB backgrounds in modes at low 

• BB backgrounds in modes at high 

• qq backgrounds in e+e“ modes at low 

• qq backgrounds in e+e“ modes at high g^; 

• qq backgrounds in modes at low g^; 

• qq backgrounds in modes at high g^, 

where low (high) g^ is defined as the mass-squared region 
below (above) the vetoed J/ijj region. In order to treat 
the control sample equivalently to the e+e” and 

datasets, we similarly train four BDTs for BB and 
gg background suppression in the low and high g^ regions, 
using a high-statistics sample of simulated B —> 
events. The pi~ BDTs are used to characterize the 
K*h^pE dataset. 

Each of the above BDTs uses a subset of the following 
observables as its input parameters: 

• the B candidate AE\ 


• the ratio of Fox-Wolfram moments i ?2 [13 

the ratio of the second-to-zeroth angular moments 
of the energy flow L 2 /T 0 [13]) both of which are 
event shape parameters calculated using charged 
and neutral particles in the GM frame; 

• the mass and AE of the other B meson in the event 
computed in the laboratory frame by summing the 
momenta and energies of all charged particles and 
photons that are not used to reconstruct the signal 
candidate; 

• the magnitude of the total transverse momentum 
of the event; 

• the probability of the vertex fitted from all the 
B candidate tracks; 

• the cosines of four angles, all defined in the GM 
frame: the angle between the B candidate momen¬ 
tum and the beam axis, the angle between the event 
thrust axis and the beam axis, the angle between 
the thrust axis of the rest of the event and the beam 
axis, and the angle between the event thrust axis 
and the thrust axis of the rest of the event. The 
thrust T of an event comprised of N particles, or 
analogousl y fo r a subset of particles in an event, is 
defined as [13 

N 

E \Pi • 

rji 2 — 1 

“ N ’ 

E \A\ 

i=l 

where the thrust axis t maximizes the magnitude 
of the thrust T, up to a two-fold ambiguity in di¬ 
rection (forward and backward are equivalent). 

As an example. Fig. |3] shows histograms of BDT output 
normalized to unit area for simulated KgTT^e'^e~ and 
K°n+ E^ p signal and combinatorial background events 
in the q\ bin. The BDT outputs for the other final states 
and g^ bins demonstrate similar discriminating power. 

Backgrounds from B —> D(—hadronic de¬ 
cays occur if two hadrons are misidentified as leptons, 
which happens at a non-negligible rate only in dimuon 
final states. These events are vetoed by requiring the in¬ 
variant mass of the K*Tr system to be outside the range 
1.84 — 1.90 GeV/c^ after assigning the pion mass hypoth¬ 
esis to the muon candidates. Residual muon misiden- 
tification backgrounds remaining after this selection are 
characterized using the K*h^pE dataset. 

For the last steps in the event selection, we adopt 
(a) the AE regions used in our recent related analy¬ 
ses of rates and rate asymmetries in exclusive B —> 
and inclusive B Xs£^£~ decays [lH l36l|- 
—0.1(—0.05) < AE < 0.05 GeV for e+e” modes; 

and (b) the gg BDT > 0.4 selection used in the inclusive 
B —>■ Xs£~^£~ analysis [13 . After all other selection cri¬ 
teria have been imposed, this gg BDT selection removes 





0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 


BDT Output 

(a) e'*‘e~ BDT output for BB background suppression in 
B+ ->• K^TT+e+e-. 



(c) BDT output for BB background suppression in 

B+ K^it+ 11+ fi-. 


FIG. 3: BDT outputs normalized to unit area for simulated 



(b) e+e BDT output for qq background suppression in 
B+ -t K'^n+e+e-. 



(d) fi+fi BDT output for qq background suppression in 
B+ ->■ K^tt+ fi+ II-. 

(solid blue line) and background (red dashed line) qf events. 


most qq background events with only trivial decreases in 
signal efficiencies. 

At the conclusion of the event selection process, some 
events have multiple reconstructed B candidates which 
typically differ by one charged or neutral pion in the 
hadronic system. The signal candidate multiplicity aver¬ 
aged across final states and bins is ~ 1.4 (^ 1.1) candi¬ 
dates per event in dielectron (dimuon) modes. In events 
with multiple signal candidates, the candidate with the 
AE value closest to zero is selected. 


III. ANGULAR OBSERVABLES EXTRACTION 
METHOD 

A. General Strategy 

We extract the angular observables and Afb from 
the data using a series of likelihood (LH) fits which pro¬ 
ceed in several steps: 

• In each q'^ bin, for each of the five signal modes 
separately and using the full tties > 5.2 GeV/c^ 


dataset, an initial unbinned maximum LH fit of 
ibeSj rn^Kn) and a likelihood ratio (Lk, defined 
below in Eq. IS that discriminates against random 
combinatorial BB backgrounds is performed. Af¬ 
ter this first fit, all normalizations and the toes- 
dependent, TO(A'7r)-dependent and L/j-dependent 
probability density function (pdf) shapes are fixed. 

• Second, in each bin and for each of the hve signal 
modes separately, ttieSj m{K'K) and Lb pdfs and 
normalizations are defined for mss > 5.27 GeV/c^ 
events (the “wes angular fit region”) using the re¬ 
sults of the prior three-dimensional fits. Only toes 
angular fit region events and pdfs are subsequently 
used in the fits for Fl and Afb- 

• Next, cos 9k is added as a fourth dimension to the 
likelihood function, in addition to toes, rn{KTr) and 
Lb, and four-dimensional likelihoods with Jfr, as 
the only free parameter are defined for toes angular 
fit region events. As above, each q^ bin and each 
of the five signal modes has its own separate 4-d 
LH function. However, a common value of Fb is 
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shared among all of the 4-d LH functions in any 
given bin. Thus, by combining LH functions 
from multiple final states, it becomes possible to 
extract and Afb for arbitrary combinations of 
the five final states here. In particular, we quote 
results using three different sets of our five signal 
modes: 

— —>■ , comprised of 

B+ K^TT+)fj.+fi-, 

B+ i4r+7r°)e+e“, 

B+ ^ K*+{-^ K°Tr+)e+e-, 

— B^ ^ , comprised of 

B° K+Tr-)n+n-, 

B° K+Tr-)e+e-. 

- B ^ K*£^£~, comprised of 

B+ K°TT+)fj.+fj,-, 

B° K+Tr-)n+n-, 

B+ K+Tr°)e+e~, 

H+ i4r*+(-> i4:°7r+)e+e“, 

B° K+Tr-)e+e-. 

• In the final step, we use the fitted value of F;, from 
the previous fit step as input to a similar 4-d fit for 
Afb, in which cos 0i replaces cos 9k as the fourth 
dimension in the LH function, in addition to toeS; 
m{KTr) and L^- 

As mentioned above, we define a likelihood ratio 
as the third dimension in the initial fit. 


TABLE II: Final signal efficiencies in the tubs angular fit 
region by mode and bin. 


Lr = 


V, 


Slg 


Slg 


Fbkg 


(3) 


where Fsig and Fbkg are probabilities calculated from the 
BB BDT output for signal and BB backgrounds, respec¬ 
tively. Fsig and Fbkg are modeled using several differ¬ 
ent functional forms depending on bin and final state. 
After the multiple candidate selection described at the 
conclusion of the preceding section and before fitting a 
dataset, a final requirement of Lr > 0.6 is made. This 
drastically reduces the number of background events at 
the cost of a relatively small loss, dependent on final state 
and bin, in signal efficiency. Table [III shows final sig¬ 
nal efficiencies in the tties angular fit region for each final 
state and q^ bin. 

The initial 3-d fit is an unbinned maximum likelihood 
fit with minimization performed by MINUIT [s^- Each 
angular result is subsequently determined by direct con¬ 
struction and examination of the negative log-likelihood 
(NLL) curves resulting from a scan across the entire Fr 
or Afb parameter space, including unphysical regions 
which provide a statistically consistent description of the 
data. 


Mode 

9o 

Qi 

ql 

93 

94 

95 

AsTt+p+p" 

0.143 

0.130 

0.146 

0.145 

0.143 

0.108 

rA 

0.184 

0.152 

0.185 

0.195 

0.194 

0.157 

K^'K^e'^ e~ 

0.121 

0.105 

0.124 

0.121 

0.110 

0.075 


0.182 

0.160 

0.185 

0.174 

0.151 

0.109 


0.230 

0.195 

0.233 

0.229 

0.209 

0.151 


B. Event Classes 

We characterize tues, m(K7r), Lr, cos 9k and cos 9i 
probability density functions in our likelihood fit model 
for several classes of events: 

• correctly reconstructed (“true”) signal events; 

• misreconstructed (“crossfeed”) signal events, from 
both the five signal modes as well as from other 
b —>• s£~^£~ decays; 

• random combinatorial backgrounds; 

• backgrounds from J/'0 and decays which es¬ 

cape the dilepton mass veto windows; 

• for the fJ-~ modes only, backgrounds from 
hadronic decays in which there is muon misidenti- 
hcation of hadrons (this background is negligible in 
e“'"e“ final states due to the much smaller, relative 
to muons, electron misidentification probability). 


1. True and Crossfeed Signal Events 

True signal events have all final state daughter particles 
correctly reconstructed. The true signal normalization 
for each final state in each q^ bin is a free parameter in 
the initial 3-d fits. For each final state, the rriES signal pdf 
is parameterized as a Gaussian with a mean and width 
fixed to values obtained from a fit to the vetoed J/ip data 
events in the same final state. Similarly, for the resonant 
K* lineshape in each final state, the signal m{KTr) pdf 
uses a relativistic Breit-Wigner (BW) with width and 
pole mass fixed from the vetoed J/if data events in the 
same final state. True signal Lr pdfs for each final state 
in each q^ bin are derived from simulated signal events. 
There is good agreement between the Lr shapes derived 
from simulated events and the Lr shapes observed in the 
charmonium control sample data. 

Equations [1] and m showing the dependence of Fr and 
Afb on cos 9k and cos 9i respectively, are purely theoret¬ 
ical expressions which must be modified to take into ac¬ 
count the experimental acceptance. We characterize the 
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angular acceptance using simulated signal events to ob¬ 
tain parameterizations of the cos 9k and cos 9i efficiency 
for each final state in each bin. 

Signal crossfeed typically occurs when a low-energy 
7 r=*= or 7r° is swapped, added or removed from the set of 
daughter particles used to reconstruct an otherwise cor¬ 
rectly reconstructed signal candidate. There can be self¬ 
crossfeed within one signal mode, feed-across between 
two different signal modes with the same final state par¬ 
ticle multiplicity, or (up) down crossfeed from (lower) 
higher multiplicity st^ modes. Simulated signal events 
are used to model these types of decays, with normaliza¬ 
tion relative to the fitted true signal yield. Averaged over 
the five signal modes and disjoint bins — qf, the frac¬ 
tion of crossfeed events relative to correctly reconstructed 
signal decays is ^ 0.4 for events in the rriES > 5.27 GeV/c^ 
angular fit region. Generator-level variations in the pro¬ 
duction of cross-feed events are considered as part of the 
study of systematic uncertainties related to the modeling 
of signal decays. 


2. Combinatorial Backgrounds 

The largest source of background is from semileptonic 
B and D decays, where leptons from two such decays 
and a K* candidate combine to form a B candidate. 
The rriES pdf for the combinatorial background is mod¬ 
eled with a kinematic threshold function |38l | whose single 
shape parameter is a free parameter in the fits. Events in 
the lepton-flavor violating (LEV) modes which 

are forbidden in the SM and for which stringent exper¬ 
imental limits exist [l^, are reconstructed and selected 
analogously to the final event selection in order to char¬ 
acterize the combinatorial background m{KTr) and Lk 
pdfs. We obtain the angular pdfs for the combinatorial 
backgrounds in the tties angular fit region using events in 
the Toes sideband region 5.2 < tues <5.27 GeV/c^. The 
LEV events additionally provide an alternative model for 
the combinatorial angular pdfs, which is used in the char¬ 
acterization of systematic uncertainties in the angular 
fits. 


3. Charmonium and Other Physics Backgrounds 

Some misreconstructed charmonium events escape the 
charmonium vetoes and appear in our q^ bins. This typ¬ 
ically occurs through bremsstrahlung by electrons, fol¬ 
lowed by incorrect recovery of the missing energy. The 
pdfs for this residual charmonium background are mod¬ 
eled using simulated charmonium signal events. 

In order to use the vetoed charmonium events as a data 
control sample, we construct a set of pdfs equivalent to 
those used in the B angular fits but which 

are appropriate for J/ifr and 1^(23) events inside, rather 
than outside, their respective vetoed mass windows. The 


BDTs in the low (high) q"^ bin are used to calculate Lk 
for events within the J/'0 {'ip{2S)) mass window. 

Gamma conversions from B K*^ events and Dalitz 
decays {n°,r]) —>■ e+e “7 of hadronic B decay daughters 
give rise to small backgrounds in However, since less 
than a single event from these sources is expected in the 
final angular fits, we do not include them in our fit model. 


4- Muon Misidentification Backgrounds 

In dimuon modes only, some events pass the final selec¬ 
tion but have misidentified hadron(s) taking the place of 
one or both muon candidates. To model these events, we 
follow a procedure similar to that described in Ref. [s^ 
by selecting a sample of events requiring that 

the candidate be identihed as a muon and the 
candidate fail identification as an electron. Using weights 
obtained from data control samples where a charged par¬ 
ticle’s species can be identified with high precision and 
accuracy without using particle identification informa¬ 
tion, the dataset is weighted event-by-event to 

characterize expected contributions in our hts due to the 
presence of misidentified muon candidates. The pdfs for 
these events are implemented as a sum of weighted his¬ 
tograms, with normalizations obtained by construction 
directly from the weighted control sample data. 


C. Initial m-ES, m{K'n) and Ln Fit 

As discussed above, the initial three-dimensional fits 
to mES, m^Kn) and L/j are done using events in the 
full Toes > 5.2 GeV/c^ range; each final state in each q^ 
bin is separately fit in order to establish the normaliza¬ 
tions and pdf shapes subsequently used in extracting the 
angular observables from the tties >5.27 GeV/c^ angu¬ 
lar fit region. Table |TTT] gives the resulting fitted signal 
yields along with statistical uncertainties for the three 
different combinations of particular final states for which 
the angular observables are extracted. As examples of 
typical fits, Eig. |4] shows fit projections in each of the 
three initial fit dimensions for A'+Tr^e+e” and 

—>■ Ar+7r“/i+p,“ in the q'^ bin. Validation of the ini¬ 
tial 3-d fit model is done using events in the J/'0 and 
'tjj{2S) dilepton mass veto windows, where we find good 
agreement between our fit results and the nominal PDG 
values for the B J/ipK* and B —>■ il;{2S)K* branching 
fractions [2^ into our final states. 


D. Angular Fit Results 

Prior to fitting the B —)> K*t^£~ angular data, we 
validate our angular ht model by using it to extract the 
K* longitudinal polarization Fl for B —>■ J/ipK* and 
B ■ijj{2S)K* decays into our signal final states, and 
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TABLE III: Fitted signal yields with statistical uncertainties. 


Mode 

lo 

gi 

ql 

<73 

ql 

Is 

B -A K*e+r 

B+ K*+£+£- 
B° -A K*°£+£- 

40.8 ±8.4 
17.7 ±5.2 
23.1 ±6.6 

31.7 ±7.1 

8.7 ±4.1 
22.9 ±5.8 

11.9 ±5.5 

3.8 ±4.0 

8.1 ±3.8 

21.3 ±8.5 

7.7 ±5.6 

13.7 ±6.4 

31.9 ±9.2 

9.0 ±4.8 
22.8 ±7.8 

33.2 ± 7.8 
9.4 ±4.2 

23.8 ±6.6 



FIG. 4: Initial 3-d fit projections for —>■ K^'K~e^e~ (top row) and ■k ~(bottom row) in q^. The plots show 

the stacked contributions from each event class: combinatorial (magenta long dash), charmonium (black dots), crossfeed (red 
short dash), total pdf (solid blue) and, in the bottom row of plots only, muon mis-identification (blue dash dots). The signal 
pdf is represented by the area between the dash red and solid blue lines. 


comparing our results to previously reported PDG val¬ 
ues [^. We also perform similar validation fits ior Afb^ 
which is expected in the SM to approach zero for lepton 
pairs from B decays to final states including charmonia. 
Recalculating the PDG averages after removing all con¬ 
tributing BaBAR results, we find no significant deviations 
from the expected values in any individual final state or 
for the particular combinations of final states used in our 
main analysis. 

Having validated our fit model with the vetoed charmo¬ 
nium events, we proceed to the extraction of the angular 
observables in each bin. Our results are tabulated in 
Tables HVl and IVl Figs. [S] and [S] show the B'^ -A 
and B° -A cos 9 k and cos0^ fit projections in Qq 

and gf. Fig.[7]graphically shows our Fl and Afb results 
in disjoint bins alongside other published results and 
the SM theory expectations, the latter of which typically 
have 5-10% theory uncertainties (absolute) in the regions 


below and above the charmonium resonances. Fig.|8]sim- 
ilarly compares the qQ results obtained here with those 
of other experiments and the SM theory expectation. 

E. Systematic Uncertainties 

We describe below the systematic uncertainties in the 
angular results arising from: 

• the purely statistical uncertainties in the parame¬ 
ters obtained from the initial 3-d mES,™(R"7r) fit 
which are used in the angular fits; 

• the Ff statistical uncertainty, which is propagated 
into the Afb fit; and 

• the modeling of the random combinatorial back¬ 
ground pdfs and the signal angular efficiencies. 
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TABLE IV: Fl angular fit results with, respectively, statistical and systematic uncertainties. 



B+ K*+£+£- 

B° K*°£+£- 

B K*£+£- 

90 

1 n r\r+0.09+0.02 

-i-u.uo_o.io_o.io 

j_n A 

-|-U.^O_0.i3_0.02 

_i_n OA 

^U.^^_0.o8-0.02 

Qi 

^ rjQ+0.18+0.09 
'^•'^^-0.13-0.14 

I r, oa+0.15+0.15 

^U.O^_0.io_o.02 

_i_n oQ+0-09+0.13 
^U.zy_0.i2_0.05 

92 

r\ Q/|+0.27+0.18 

u.zq:_o.39_o.io 

1 fl 1 o+0.16+0.02 

^U.iO_o.i2_0.10 

in 1 '7+0. 14 +0.02 
^U.i/_0.15_0.02 

93 

1 n ir+0.14+0.05 
i-U.10_o.i3_o. 08 

+U.4»_o.i6_o.05 

^U.OU_0.11 _0.07 

94 

-LO nt^i+o-^^+o.ie 
-hU.UO_o.i6_o.15 

in ac;+oo9+o.06 
^U.^O_o.i4_o.06 

in oa+0.15+0.07 

^U.O^_0.10_0.10 

95 

—LO '70+0.20+0.10 
-|-u. /Z_o 31_0.21 

+(J.4t5_0.i2_0.11 

1 n t::o+0.10+0.07 
^U.OO_o.i2_0.14 


TABLE V: Afb angular fit results with, respectively, statistical and systematic uncertainties. 



B+ K*+£+£- 

B° K*°£+£- 

B K*£+£- 

9o 

_lA qo+0.18+0.08 

-|-U.OZ_0.i8_0.05 

j_n nf\+o+^+o o6 

-hU.UO_o.i8_o.05 

_i_n 91 +0.10+0.07 
^U.zi_0.i5_o.09 

9i 

_L0 44+0-20+0.13 

-|-U.^O:_o 22-0.16 

n 1 q+O- 23+0.10 
— U.iZ_o 21 -0.21 

_i_n in+o.i6+o.o8 
^U.iU_o.i5_o.i9 

92 

1 n '7n+0-21+0.36 
-|-U. /U_o.38_o.49 

in oq+0.21+0.12 

^U.OO_o.30-0.11 

in 44+0.15+0.14 

+U.44_0.i8-0.11 

93 

1 n 11+0.22+0.08 
-|-U.11 _o.28-0.20 

in 1'7+0.14+0.08 

^U.i/_o.i6-0.08 

in 1 c:+0.14+0.08 
^U.iO_o.i2-0.05 

94 

1 n O1+0.32+0.11 
-rU.zi_o 33_o.24 

in 40+0.12+0.17 

^U.^U_o.i8_o.16 

-1-0 +0.14 

^U.^Z_o i7_o.i3 

95 

1 n An+o-26+0.18 
-|-U.^U_o.21_o.17 

in 9Q+0.14+0.10 
^U.zy_0.i7_o.io 

_i_n 9Q+0.07+0.10 
^U.Zi:;_ 0 .io_ 0.12 


We additionally examined several other possible 
sources of systematic uncertainty, but found no signifi¬ 
cant contributions due to: 

• modeling of the signal crossfeed contributions to 
the angular fits; 

• the parameterization of the signal Gaussian toes 
and resonant m{KTr) shapes that are extracted 
from the relatively high-statistics J/'i/' control sam¬ 
ples; 

• possible fit biases which, to relatively very good 
precision, were not observed in any of the data con¬ 
trol sample angular fits; 

• characterization of toes peaking backgrounds from 
muon mis-identification and charmonium leakage; 

• variations in event selection. 

We combine in quadrature the individual systematic 
uncertainties to obtain the total systematic uncertainty 
on each of the angular observables; these are are given 
in Table s which is placed after the detailed discussion 
below for each family of systematic uncertainties. 

In the initial fits that determine the signal yields, we 
allow the random combinatorial toes shape and normal¬ 
isation, as well as the signal yield, to float. We then 


fix these parameters at their central values for the an¬ 
gular fits. To study the systematic uncertainty associ¬ 
ated with these fixed parameters, we vary each param¬ 
eter from its central value by its ±lcr statistical uncer¬ 
tainty, accounting for correlations among the fit param¬ 
eters, and then redo the angular fit. To control for sys¬ 
tematic fit results that deviate from the nominal central 
value mainly from statistical effects rather than system¬ 
atic ones, we additionally examine fit results obtained 
from ±(0.8,0.9,1.1, 1.2 )ct variations. These small varia¬ 
tions on the ±1(7 values should also result in similarly 
small variations, in the absence of any statistical effects, 
on a ±1(7 systematic fit result. For the bulk of the sys- 
tematics, where the series of fit results for each of the ad¬ 
ditional variations is linearly distributed around the mid¬ 
dle 1(7 fit result, the la variation is considered robust. In 
the relatively few cases where the disagreement between 
the nominal la variation and the value of the la variation 
interpolated from the additional (0.8, 0.9,1.1,1.2)(7 varia¬ 
tions is statistically significant, the interpolated la value 
is used to assign the systematic. All deviations from the 
nominal fit central value are then added in quadrature to 
obtain the overall systematic uncertainty attributable to 
this source, which is given in Table IVII 

The cos 9 k fit yields the central value and statistical 
uncertainty for Fl in each bin, which is subsequently 
used in the fit to the cosOi distributions to extract Afb- 
To study the systematic uncertainty on Afb due to the 
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(a) cos6»x (b) cos0^ gg 




FIG. 5: —>■ K*'^£'^£ angular fit projections. The shaded areas show the contribution to the total fit from each individual 

final state: (green vertical lines) —>■ (red diagonal lines) B'^ —>■ K^'R^e^e~\ (blue horizontal lines) B'^ —>■ 

K+Tv^e+e- '. The overlaid dashed line shows the total signal contribution summed over the three individual final states. Each 
colored band includes both signal and background events in a given final state. 


purely statistical uncertainty, we vary the value of 
Fl by ±lcr from its fitted value, and redo the cos9i fits 
with the new value of We determine the systematic 
uncertainty from the shift in the central value of Afb rel¬ 
ative to the nominal fit for ±lcr variations of these 
are given in Table IVTIl As with the variations described 
in the preceding paragraph, additional fits for several F^ 
variations surrounding the nominal ±lcr values are per¬ 
formed. We then apply the same quality criterion as for 
the preceding systematic and, where this criterion is not 
met, assign the Fl systematic using an interpolated Ict 
value rather than the fitted Ict variation. 

The angular combinatorial background shapes are 
derived from the toes sideband region and are non- 
parametrically modeled directly from these data. We ex¬ 
amine several variations on the modeling, and addition¬ 
ally use the LFV events (described above) as an alterna¬ 


tive dataset from which the angular background pdfs are 
drawn. We assign a systematic uncertainty associated 
with the modeling of these pdfs by using 20 different 
variations of the non-parametric modeling and refitting 
for Fl and Afb- We take the largest of the deviations 
between the default nominal fit and these varied fit re¬ 
sults, and to this add in quadrature the deviation from 
the nominal fit obtained using the LFV dataset; the re¬ 
sulting systematic uncertainty is given in Table IVIIII 

Finally, to study a possible systematic uncertainty on 
Ff and Afb as a function of their true physical values, 
we generated and reconstructed simulated events with 
varied values of the underlying Wilson coefficients Cr,Cg, 
and Cio in order to produce a range of near-maximal, but 
physically allowed, asymmetries. These datasets are used 
to produce signal efficiency histograms differing from the 
default ones, which use the expected SM values for the 
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FIG. 6: —>■ K*^t^l~ angular fit projections. The shaded areas show the contribntion to the total fit from each individual 

final state: (red diagonal lines) 5° K'^■, (blue horizontal lines) . The overlaid dashed line 

shows the total signal contribntion summed over the two individual final states. Each colored band includes both signal and 
backgronnd events in a given final state. 


Wilson coefficients. This allows different regions of the 
angular distributions to contribute with different rela¬ 
tive weight depending on the magnitude and sign of the 
underlying angular asymmetries. Applying these alter¬ 
native signal efficiency histograms, we measure the shifts 
in the fitted values of Fl and Afb, and assign as the 
systematic the sum-in-quadrature of each deviation from 
the nominal central value; the resulting systematic un¬ 
certainty is shown in Table llXl 


F. Extraction of P2 from the Angular Fit Results 

As mentioned above in the Introduction, Fl and Afb 
can be used to parameterize an additional angular ob¬ 
servable, P 2 = (—2/3) *Afb/(1 — Fl), which has dimin¬ 
ished theory uncertainty and greater sensitivity to non- 


SM contributions than either Fl or Afb alone [ 2 ^ [ 2 ^ . 
Table |XT] gives our results for P 2 . The 68% confidence 
intervals quoted are frequentist and derived from ensem¬ 
bles of fits to simulated datasets randomly drawn from 
the correlated confidence-level contours for Fl and Afb- 
For the disjoint mass-squared bins q\ to g|. Fig. |9]graph- 
ically shows our results overlaid on the SM expectations 
from theory (as given in Table 3 (KMPW) of Ref. Hi) 
in the mass-squared region below the J/ip. In the 
mass-squared bin, the SM expectation (from the same 
source) for P 2 is 0.11 ± 0.10, in slight tension with our 
experimental result. 
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* 

Belle 

• 

CDF 

□ 

LHCb 

0 

CMS 

ir 

ATLAS 

m 

Babar K* 

▼ 

Babar K*° 

A 

Babar K** 



(a) Fl. 



(b) Afb- 


FIG. 7: Fl (top) and Afb (bottom) results in disjoint bins, along with those of other experiments and the SM expectations 
(blue dashed lines, which also define the extent of each individual bin): (black filled star) Belle [l^. (black filled circle) 
CDF [ 2 ^, (black open square) LHCb [ 2 ^, (black open circle) CMS [ 2 ^, (black open star) ATLAS [ 2 ^, (blue filled square) BABAR 
B —>• K*£'^£~ , (red filled down-pointing triangle) B^ —>■ K*°£'^£~ , (magenta filled up-pointing triangle) B+ —>■ K*^£'^£~ . The 
BABAR 55 results are drawn in the 14 < < 16 GeV^/region, however, they are valid for the entire q^ ^ 14 GeV^/region. 




(a) Ff- (b) Afb- 


FIG. 8: On Fl (left) and Afb (right) results, along with those of other experiments and the SM expectation (vertical 

lines) [J-la, [3]- 
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TABLE VI: Angular observable systematic uncertainties from the initial 3-d fit. 



B+ K*+e+£- 

Fl systematic 

B° K*°e+i- 

B K*£+r 

B+ K*+£+£- 

Afb systematic 
B° K*°£+£- 

B -> K*£+£- 

90 

+0.02 -0.09 

+0.02 -0.02 

+0.02 -0.02 

+0.05 -0.04 

+0.01 -0.04 

+0.02 -0.07 

9i 

+0.09 -0.13 

+0.02 -0.02 

+0.02 -0.05 

+0.12 -0.08 

+0.07 -0.02 

+0.07 -0.08 

92 

+0.18 -0.05 

+0.02 -0.01 

+0.02 -0.02 

+0.34 -0.48 

-0.02 -0.08 

+0.09 -0.07 

93 

+0.05 -0.07 

+0.02 -0.02 

+0.05 -0.06 

+0.02 -0.19 

-0.02 -0.04 

+0.01 -0.02 

94 

+0.11 -0.14 

+0.02 -0.06 

+0.02 -0.10 

+0.09 -0.23 

+0.15 -0.11 

+0.13 -0.10 

95 

+0.02 -0.19 

+0.02 -0.10 

+0.02 -0.14 

+0.16 -0.09 

+0.05 -0.02 

+0.08 -0.02 


TABLE VII: Systematic uncertainty in Afb from the experimental determination of Fl- 



B+ -> K*+£+£ 

Afb systematic 
- ^ K*°£+£- 

B K*£+£- 

9o 

±0.04 

±0.04 

±0.04 

9i 

±0.04 

±0.07 

±0.04 

92 

±0.07 

±0.07 

±0.08 

93 

±0.03 

±0.06 

±0.04 

94 

±0.04 

±0.07 

±0.06 

95 

±0.08 

±0.07 

±0.07 


TABLE VIII: Systematic uncertainties from combinatorial background modeling. “—” denotes where there is no uncertainty 
associated with a particular systematic. 



B+ K*+e+£- 

Fl systematic 

K*°£+£- B K*£+£- 

B+ K*+£+£- 

Afb systematic 

B° K*°e+£- 

B K*£+£- 

9o 

— -0.05 

— — — — 

+0.04 — 

— — 

— -0.04 

9i 

±0.02 -0.02 

— — — — 

+0.05 — 

— — 

— — 

92 

— -0.05 

— — — — 

— -0.07 

— -0.04 

— — 

93 

— — 

— — — — 

— — 

— — 

— — 

94 

+0.10 — 

— — — — 

— -0.04 

— — 

— — 

95 

— -0.10 

— -0.05 — — 

+0.04 -0.08 

+0.04 — 

— — 
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TABLE IX: Systematic uncertainties from signal angular efficiency modeling. “—” denotes where there is no uncertainty 
associated with a particular systematic. 




Fl systematic 


Afb systematic 



B+ ^ k*+£+£- 

B° K*°£+£- B K*£+£- 

B+ k*+£+£- 

B° K*°£+£- 

■ B K*£+£- 

1 

Oto 

— - 0.02 

— — — — 

+0.04 — 

+0.04 — 

+0.05 — 

9i 

+0.02 -0.04 

+0.14 — +0.13 — 

— -0.13 

— - 0.20 

— -0.17 

92 

— -0.07 

— - 0.10 — - 0.02 

+ 0.12 — 

+0.09 — 

+0.07 — 

93 

— -0.04 

+0.04 -0.05 +0.02 -0.04 

+0.08 — 

+0.06 -0.04 

+0.07 -0.02 

94 

+0.07 -0.05 

+0.06 — +0.07 — 

+0.06 — 

— -0.09 

+ 0.02 -0.06 

95 

+ 0.10 — 

+0.02 — +0.07 — 

— -0.09 

— -0.08 

— - 0.10 

TABLE X: Total systematic uncertainties. 



Fl systematic 


Afb systematic 



B+ K*+e+£- 

B° K*°£+£- B K*£+£- 

B+ -F K*+£+£- 

B° K*°£+£- 

■ B K*£+£- 

1 

Oto 

+ 0.02 - 0.10 

+ 0.02 - 0.02 + 0.02 - 0.02 

+0.08 -0.05 

+0.06 -0.05 

+0.07 -0.09 

9i 

+0.09 -0.14 

+0.15 -0.02 +0.13 -0.05 

+0.13 -0.16 

+ 0.10 - 0.21 

+0.08 -0.19 

92 

+0.18 - 0.10 

+ 0.02 - 0.10 + 0.02 - 0.02 

+0.36 -0.49 

+ 0.12 - 0.11 

+0.14 -0.11 

93 

+0.05 -0.08 

+0.05 -0.05 +0.05 -0.07 

+0.08 - 0.20 

+0.08 -0.08 

+0.08 -0.05 

94 

+0.16 -0.15 

+0.06 -0.06 +0.07 -0.10 

+0.11 -0.24 

+0.17 -0.16 

+0.14 -0.13 

95 

+ 0.10 - 0.21 

+0.02 -0.11 +0.07 -0.14 

+0.18 -0.17 

+ 0.10 - 0.10 

+ 0.10 - 0.12 



TABLE XI: P 2 results with total uncertainties. 




B+ ^ K*+£+£- B 

° ^ K*°£+£- 

B K*£+t 



9o 

- 0 . 22 «;lt 


-0.18t°;« 



9i 

-0.29l°;l? 

+0.12t°S, 

-0.09t°;« 



92 

-0.38l°:i 

A 97+0.25 
'^•^'-0.24 

-0.35t°;j^ 



93 

-0.09^0.2^ 

-0.22t°of2 




94 

—0 1 

-0A8+°f, 

— 0 49+0-26 

U.+^_ 0.20 



95 

U.yO_0.96 

A 07 + 0.28 
-0.24 

— 0 41 +0-34 
^•+-*-- 0.21 



IV. CONCLUSION 

the charmonium 

1 resonances. Fig. | 8 ] similarly compares 




the Qq results obtained here with those of other exper¬ 
iments and the SM theory expectation. As shown in 

In conclusion, we have measured in bins of dilepton 

these figures, our B^ —^ K*'^£~^£~ results are in rea- 

mass-squared the fraction F^ of longitudinally polarized 

sonable agreement with both SM theory expectations 

K* decays and the lepton forward-backward asymmetry 

and other experimental results. 

Similarly, although with 

Afb in the decays 5+ —)■ 

K*+£+£-, B° K*°£+£- 

relatively larger 

uncertainties, - 

we observe broad agree- 

and B —^ K*£'^£ . Results for the —)> K*'^£'^£ fi- 

ment of the B~^ K*~^£~^£ 

results with those for 

nal state are presented for the first time here. Fig. [7] 

B° K*°£+£- 

. However, in 

the low dilepton mass- 


graphically shows our Fl and Afb results in disjoint 
bins alongside other published results and the SM the¬ 
ory expectations, the latter of which typically have 5- 
10 % theory uncertainties in the regions below and above 


squared region, we observe relatively very small values for 
Fl in —>■ exhibiting tension with both the 

^ K*^£'^£~ results as well as the SM expectations. 
Similarly, as shown in Fig. in the same mass-squared 
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FIG. 9: P 2 results with total uncertainties.The blue boxes show the SM theory expectation in the low mass-squared region; 
there are no comparable calculations in the high mass-squared region, where the black lines simply denote the extent of the gl 
and gl bins. 


region, there appears to be some tension between the ex¬ 
perimental results and the expected SM values for P 2 . 
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